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Abstract. - The complexation of oppositely charged colloidal objects is considered in this 
paper as a thermodynamic micellization process where each kind of object needs the others to 
micellize. This requirement gives rise to quantitatively different behaviors than the so-called 
mixed-micellization where each specie can micellize separately. A simple model of the grand 
potential for micelles is proposed to corroborate the predictions of this general approach. 



Complexation of oppositely charged colloidal objects has been the subject of continuous 
and important researches since the seventies because of the numerous applications ranging 
from flocculants and coatings to systems mimicking biological phenomena . The interest in 
this last field has grown in recent years as the drug carrier potential of polymer complexes has 
been realized ^ : most of the non- viral vectors used for gene delivery are made of complexes 
between negatively charged DNA and oppositely charged objects that can be either surfactants 
(lipoplexes) [3] or homo or co-polymers (polyplexes) . Since a drug delivery vector has to 
be able to cross several biological and physical obstacles, the tailoring of a multifunction carrier 
challenges nowadays biologists as well as physico-chemists |S] . The common feature of these 
complexes is their reliance on electrostatic self-assembly. Complementary to the experimental 
approaches looking for various drug delivery systems based on electrostatic complexation, 
there is obviously a need for more thorough understanding of the underlying mechanisms of 
self-assembly and of the relevant parameters for controlling them . 

In this Letter I present the general formalism required to describe the thermodynamic 
properties of micelles formed by complexation of oppositely charged objects and having a 
well defined gaussian size distribution. This allows us to address in particular the role of 
the mixture stoichiometry on the micelle compositions. In order to derive the most general 
features of the complexation viewed as a micellization process, the kind of objects considered 
will not be specified . The word "micellization" is used here in a generalized way: monodis- 
perse amphiphilic molecules self-assemble forming micelles above the so-called critical micelle 
concentration as they are not soluble as a single specie 0. In the case of complexation, both 
components of the complexes are water-soluble, but they still form aggregates that have been 
shown in the case of various polymer systems to exhibit a core-corona structure similar to the 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



polymeric micelles made of neutral hydrophobic-hydrophilic copolymers 2 . This extension of 
the micelle concept from aggregates of amphiphilic molecules to electrostatic complexes has 
been put forward quite recently by Cohen-Stuart et al. I will emphasize in this work the 
major differences between these two kinds of micellization 9^. As an example, a phenomeno- 
logical model of grand potential describing micelles made of oppositely charged objects is 
proposed and the phase behavior of such a solution is derived numerically. 

Depending on the kind of systems considered and on the stoichiometry of the mixture, two 
main features are observed pj: one can have either a macroscopic phase separation, associated 
with the neutralization of charges, or a partial micellization where each aggregate carries a net 
charge. For the sake of simplicity, I will focus mainly on the partial aggregation in this work . 
Physical realizations of such systems are, for example, complexes formed with polyelectrolytes 
and so-called double hydrophilic polymers - diblock copolymers containing a polyelectrolyte 
block of opposite charge and a water-soluble neutral block (typically polyethylene glycol or 
PEG) pi ll0l[TH ll2 | . The macroscopic aggregation is limited by the neutral blocks forming a 
corona around a complexed core. This choice of systems allows us to focus on the effect of 
the stoichiometry of the solution on the properties of the micelles (13) . Moreover, complexes 
of classical oppositely charged polyelectrolytes are shown to aggregate only partially for non 
neutral stoichiometrics of the mixture and therefore some results obtained here might be also 
useful in this case. For strongly charged polyelectrolytes, the complexation is thought to 
take place out of thermodynamic equilibrium. Nevertheless, the knowledge of equilibrium 
properties is the first step towards kinetic or dynamical modeling. 

Let us consider a bidisperse solution made of oppositely charged objects. Since I do not 
allow for macroscopic phase separation in the system, the small ions neutralizing the solution 
are not explicitly taken into account, to simplify the discussion; the effect of the salt can be 
thought to be incorporated into the free energy of the micelles. Each object is supposed to 
carry a ±1 charge. The energy unit is set to kT = 1 for the remaining of the paper. The free 
energy density describing the solution is given in the dilute limit by 

oo 

F ^c+[logc++F+] + epilog C-+F^]+ ^ Cp^,p_[logCp^^p^ + Fp^^p_] (1) 

p+>p-=i 

The quantities c+ , c_ , Cp_^ ^p_ are respectively the concentrations of positive and negative un- 
aggregated objects and the concentration of micelles made of p+ and p_ such objects. The 
corresponding excess chemical potentials are F^,F- and i^p+,p_. This free energy also takes 
into account the translational entropy of single objects and micelles. The equilibrium state of 
the system is determined by its minimization with respect to all concentrations, and subject 
to the following mass conservation equations 

oo 

C++ P+^P+^P- = (2) 
oo 

c- + X! P-^P+.P- = (3) 

The total concentration of positive and negative objects are respectively and (/>_. In- 
troducing the two Lagrange multipliers /i+ and /i_ corresponding to the mass conservation 
equations, the minimization equations can be written in the form: 



c+ = exp{-(F+ - 



(4) 
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c_ = exp{-(F_~M-)} (5) 
Cp+,p_ = exp{-(Fp^,p_ - - (6) 

/z+ and /i_ are interpreted as the chemical potential of positive and negative objects. The last 
equation shows how the concentration of micelles is related to the grand potential Qp^ p_ = 

Fp+,P- -P-f^-- 

In the case studied here of well defined micelles having a gaussian size distribution, the 
grand potential has a minimum located at the finite values p+ and p_ above some critial value 
for the chemical potentials. Therefore the sums in the mass conservation equations can be 
approximated by a steepest descent method. This amounts to considering an effective two 
state model for the components in the solution, the unaggregated object and the micelle state. 
The mass conservation equations can then be rewritten 

l + cxp{-np_^^p_}exp{F+ - fi+}p+A = cxp{-{^i+ - nl)} (7) 
l + exp{-r2p_|_,p_}exp{F_ -^_}p_A exp{-(/i_ - ^" )} (8) 

with the following notations A = 27r/y^(5+(5_ — Sj., Si id'^Fp^^p_/dpf)p^^p_ and d± = 

{d'^ Fp+,p- /dp+dp ) taking into account the width of the micelle size distribution. Each 
total concentration has been rewritten in the preceding equations by introducing the chemical 
potential of objects without any association /i° = log (pi + Fi. 

The solutions of EqEI and |S1 for /i+ and /i_ together with the minimization equations 
for the grand potential allows us to determine all the equilibrium properties of the solution, 
e.g. to draw its phase diagram, once the grand potential of the micelles is known. One 
can understand the micellization behavior in the solution by following particular paths in 
the phase diagram. This is done numerically below with a crude model of grand potential 
bearing the two essential features of the particular micellization studied here: one requires 
that (i) micelles cannot form if only one component is present in solution and that (ii) any 
net charge of the micelles has an energetic price. The first requirement can be easily taken 
into account by an attractive term of the form Fatt = Ap'^p^_. For positive values of A, the 
sum of the exponents has to be lower than 1, so that this form describes an attraction. This 
can be thought as a "surface tension" term favoring the growth of the complex, similar to 
the micelles made of regular diblock hydrophilic/hydrophobic copolymers, although strictly 
speaking such a term would be proportional to the two-third power of the sum of aggregation 
numbers ( which is unable to fulfill requirement (i)). Using nevertheless this analogy, the 
exponents are chosen to be a = /? = 1/3 for a symmetric system, so that this term favoring 
aggregation corresponds exactly to the single micellization case in the case of zero net charge. 
Note that the precise values of the exponents do not change the qualitative behavior described 
below. The second requirement can be taken into account with a quadratic energetic price 
for net charges, Fasym — B(j>^ — Such a form occurs for example simply by taking 

into account the loss of translational entropy by small ions neutralizing the excess charge in 
the complex |lHlll4j . Finally, the prohibition of macroscopic phase separation through the 
complexation is ensured by a term representing for example the energetic cost of a polymeric 
corona around the complexed core, Fcorona = Cip+ +P-)'^^'^ JS]. The choice of a corona-like 
stabilization of the complex is obvious for cases where there is actually a corona, while it 
might look odd for a mixture of oppositely charged homopolyelectrolytes. Nevertheless, as 
proposed in reference ^3] for diblock polyampholytes, the excess charge might be favorably 
located on arms surrounding the neutral complexed core. 

The grand potential is therefore written as 

f^P+,p_ = Ap'fp]/'' + B{p+ -p_f -f C{p+ +p-f'^ - fi+p+ - iji-p- (9) 
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Fig. 1 - (a) Chemical potential of positive (circles) and negative (crosses) objects as a function of 
the bare chemical potential of positive objects /^^j. for /I'L — 15. The dashed line represents fi+ = fj!^. 
(b) Particular paths in the phase diagram for different starting points (also indicated by arrows): 
/u^ = 11 (pluses), fj,^ = 13 (crosses), fj,^ — 15 (circles). The parameters of the grand potential have 
been set equal to A = 30, B = 10, C = 1, F+ = F_ = 0. 



The mass conservation equations Eq[7| and |S1 are then solved numerically for a given set of 
parameters together with the minimization equations of the grand potential to obtain the 
chemical potentials of positive and negative objects. As will be shown below, the results can 
be traced back to the general properties of this specific kind of micellization. The properties 
of the solution can be probed through particular paths in the phase diagram. Experimentally 
relevant paths correspond for example to the increase of the number of positive objects at 
fixed number of negative objects. An example of numerical solution for the chemical potentials 
along such a path is given on figure^ At very low concentration of positive objects, there is 
no association and /i± — Above a critical concentration of positive objects, micelles start 
to form and the chemical potential of positive objects saturates while the chemical potential 
of negative objects remains unaffected: almost every added positive object is consumed by 
forming a complex with negative objects. This situation is very similar to the monodiperse 
micellization case of amphiphilic molecules, where the chemical potential of molecules is almost 
constant above the critical micellar concentration (CMC). The negative objects play in the 
case of the present study the role of "linkers" that help the positive objects to aggregate. The 
nearly constant value of chemical potential of positive objects, which gives also the onset of 
micellization, can be found by requiring that the term associated with the micelles in Eq[7|is 
of order one: 

^, Fp+.P- - P+ log [p+A] ^i°_p- 

— : : 7 (J-Uj 

P+-1 P+-1 

Just above this onset of micellization, most of the negative objects are still in the unaggregated 
state. The micelle term in the mass conservation equation of negative objects becomes of order 
of the isolated object term for a higher value of fuP^. Using Eqd the micelle term in the mass 
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conservation equation of negative objects is rewritten just above the onset of micellization as 

exp{-fip+,p_}exp{i^_ - M-}P-A = (11) 

Since the micelles under consideration are such that they cannot handle large non-neutral 
stoichiometric ratio ^ 1, the micelles start to dominate the behavior of negative objects 

very close to the value = yP_—F-+F+, which corresponds to the equality of concentrations 
of positive and negative objects. This particular value is due to the unit charge carried by each 
object. The straightforward generalization of this result to objects of charge g+ and q- is given 
by /i" = yP_ — F- + Fj^ + log(g+/(7_). Above this second critical value, the consumption of 
added positive objects in complexes can not be as high as in the regime of excess of negative 
objects. Therefore the unaggregated objects dominate once again the behavior of positive 
objects, and their chemical potential is asymptotically given by ~ /i^. The chemical 
potential of negative objects reaches the asymptotic value 

^_^ FP.-P--F- _^I^^_P±^. (12) 

P- P- P- 

The decrease in the chemical potential of negative objects is associated with the energetic 
gain of complex formation and the increasing number of negative objects involved. All these 
features are indeed observed for the proposed grand potential through the numerical solutions. 
Different starting concentrations of negative objects are shown in figure^to end up roughly on 
the same "micellization" line of the phase diagram /i-), which corresponds approximately 
to a constant value of the grand potential minimum 

f^P+,P_ ^ F_ + logj5_A - ^I'L (13) 

In the case of bidipserse micellization where the two components can micellize separately, 
the behavior is quantitatively different for the same type of path in the phase diagram [S]: 
above the critical micellar concentration the stoichiometry of the micelles follows the stoi- 
chiomctry of the mixture, since asymmetry in the micelle composition costs essentially noth- 
ing in this case. As the concentration of one of the components is increased, the composition 
of the micelles gets close to the one of pure micelles of this component, and both chemical 
potentials saturate. Further increase of the concentration leads to an increase in the number 
of pure micelles of the added component. Moreover it has been shown that under certain 
conditions, the mixed micellization occurs in a discontinuous way, the aggregation number of 
one component jumping from to a finite value 7^ 1 in already preformed micelles of the 
other components 9 . This leads to a bimodal distribution of micelles. Such a phenomenon 
is not predicted in the case of the micellization of oppositely charged polymers. 

Of more interest from the standpoint of the drug carrier potential of complexes are quan- 
tities like the net charge and the weight of the complexes. These are shown for the grand 
potential EqlHI as a function of the mixing ratio X = 4)^/4)^ in figure El A unit mass for 
each object has been chosen for the sake of simplicity. The representation as a function of X 
is often used to plot experimental data on complexes. As shown on figure [3 the net charge 
and the weight of the complexes depend not only on the mixing ratio, but also on the overall 
concentration, a fact that is indeed observed experimentally. The largest weight is achieved 
for neutral stoichiometry and increasing overall concentration jl2l . Note that the saturation 
in the weight of the micelle, if any, is quite slow as a function of the overall concentration of 
the solution. The net charge of the micelles is shown to be tuned for the symmetric system 
studied here by the mixing ratio of the solution, although the effect is small here for the range 
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(a) (b) 



Fig. 2 - Structural parameters of the micelles as a function of the mixing ratio X — (j)+/(j>~ for 
different starting points, (a) Weight p+ + p- and (b) net charge of the micelles p+ — p-. Starting 
points: f/L = 11 (pluses), fjH. — 13 (crosses), = 15 (circles). The parameters of the grand potential 
have been set equal to A = 30, B = 10, C = 1, F+ = F_ = 0. 

of chosen parameters. For /.t*^ = 11, the system is very close to the minimal concentration 
of negative objects required to grow any complexes; this explains the slow increase of both 
the mass and the net charge. These results still show that one must be careful in modeling 
overcharging phenomena in complexes: studies like the one of Kunze and Netz on complexes 
formed by one DNA with excess flexible polycations can be naturally complemented by the 
formalism used here to address for example the effect of finite DNA concentrations 6 . 

The present study provides a rigorous framework to study the thermodynamic behavior of 
micelles made of oppositely charged polymers and in particular the effect of concentration on 
the main features of micelles. The crude choice of grand potential, based on the requirement 
that it describes in a phenomenological way species that need each other to micellize, exhibits 
the qualitative features of polyelectrolyte complexation Nevertheless, a more precise 

choice of grand potential, as well as a more thorough investigation of the relation between 
effective, i.e. A, B, C, and physical parameters (salt amount, charge fraction, etc..) is needed 
to address for structural aspects ( e.g. radii) of micelles as well as the effect of salt on the 
general properties of micelles. This is the topic of further works in progress. 
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